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Abstract 

In string/M-theory with a large compactification radius, some axion-like 
moduli can be much lighter than the gravitino. Generic moduli in gauge- 
mediated supersymmetry breaking models also have a mass far below the weak 
scale. Motivated by these, we examine the cosmological implications of light 
moduli for the mass range from the weak scale to an extremely small scale of 
0(10-26) eV, and obtain an upper bound on the initial moduli misalignment 
for both cases with and without a late entropy production. 
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I. INTRODUCTION 



A very weakly-interacting light scalar field can be cosmologically troublesome. In the 
early universe with the Hubble expansion parameter H much bigger than the scalar boson 
mass m^, this scalar field would take an initial value 50 significantly different from its vacuum 
value at present. Due to this misalignment, starts to coherently oscillate at a later time 
and subsequently its energy density behaves as a non-relativistic matter. Depending upon 
its lifetime, such a coherently oscillating scalar field may spoil the big bang nucleosynthesis 
or produce too much gamma- or X-rays through its late-time decay, or may overdose the 
universe. 

The most well-known and first studied example of such a light scalar boson is the invisible 
axion which has been introduced to solve the strong CP problem The mass and the 

typical misalignment of the invisible axion are determined by a single unknown parameter, 
the axion decay constant fa which corresponds to the scale of spontaneous f/(l)pQ -breaking, 
as 5a ^ fa and rria ~ /7r^7r//a where /,r and m,r are the pion decay constant and the 
pion mass, respectively. If there is no entropy production after the QCD phase transition, 
the requirement that the relic axion mass density does not overdose the universe leads 
to the famous constraint fa ^ 10^^ GeV. It was later realized that generic hidden sector 
supergravity model predicts a cosmologically troublesome light scalar field, the Polonyi 
field, whose typical misalignment is given by the Planck scale Mp = l/y/8nG^ ^2x 10^^ 
GeV Again in the absence of a late time entropy production after the Polonyi field 
starts to oscillate, it would completely spoil the successful big-bang nucleosynthesis unless 
the Polonyi mass ^ 40 TeV. 

This Polonyi problem has been recently revived in the context of string theory 0]. The 
(approximately) degenerate string vacua are described by moduli fields whose typical mis- 
alignments are given by either the string scale or the compactification scale which is close 
to Mp. In the previous studies, the moduli masses were assumed to be of order the grav- 
itino mass 1712/2 which would be of order of the weak scale if the supersymmetry breaking 
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is transmitted to the supersymmetric standard model sector by gravitational interactions. 
However due to approximate non-linear global U{1) symmetries in string/M-theory, some of 
the axion-like moduli in string theory can be much lighter than the gravitino In particu- 
lar, in models with a large compactification radius, the mass of such light axion-like moduli 
behaves as ^ €"'^^'^^■^^3/2 where T is the overall modulus field whose value corresponds 
to the compactification radius-squared in the heterotic string length unit. As we will discuss 
in section II, Re(T) can be as large as in the M-theory limit, implying that 

can be extremely smaller than m3/2, for instance as small as 10~^'*m3/2. 

Although not as dramatic as the axion-like moduli in string/ M-theory, generic moduli in 
models with a gauge-mediated supersymmetry breaking can also be much lighter than the 
weak scale . In gauge-mediated models, the supersymmetry breaking scale A is related to 
the weak scale Mw by some loop suppression factor: Mw = (f )'A where the integer / counts 
the number of loops involved in transmitting supersymmetry breaking to the supersymmetric 
standard model sector. Then for a reasonable value of I, the moduli mass ~ A^/Mp 
would be in the range far below the weak scale. 

Motivated by the above observations, in this paper we wish to study some cosmological 
aspects of a generic light modulus with an arbitrary mass in the range from the weak scale 
to an extremely small mass scale of 0(10^26) ^^^^^ ^-^at no meaningful cosmo- 

logical bound is obtained for the moduli mass below 10~^^ eV. In fact, moduli cosmology in 
gauge-mediated supersymmetry breaking models has been discussed recently 0, but only 
for a rather narrow mass range = 10 keV — 1 GeV. The organization of this paper goes 
as follows. In section II, we discuss in more detail the masses of the light axion-like moduli 
in string/M-theory to make our motivation more clear. In section III, we examine some 
generic features of the moduli dynamics in the early universe to see how it depends on the 
parameters involved and also identify the initial moduli misalignment 6(j). We then use the 
known cosmological and astrophysical observations to obtain an upper bound of the moduli 
misalignment 50 in section IV. In this regard, we first consider the case that there is no 
entropy production after the moduli oscillation begins and summarize the results in Figure 
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1. As is well known, the light moduli density (relative to the entropy density) can be diluted 
if there occurs an entropy production after the moduli oscillation begins. We thus finally 
discuss the dilution of the light moduli density in various possible cosmological scenarios 
with a late entropy production, including the case that the entropy-producing field is a 
massive moduli with ^ (9(40) TeV and also the case that ip corresponds to the flaton 
field triggering thermal inflation |^. The relaxed bounds on 6(f) for the cases with a late 
entropy production are depicted in Figures 2 and 3. 



II. LIGHT AXION-LIKE MODULI IN STRING/M-THEORY 

To make our motivation more clear, let us discuss in more detail the masses of axion-like 
moduli in compactified string/M-theory. As is well known, the theory predicts the dilaton 
superfield S and the Kahler moduli superfields 7} (/ = 1,2,..., /ii i) for generic compactifi- 
cations preserving the four-dimensional supersymmetry 0. The scalar components Re(S') 
and Re(Tf) of these superfields determine the four-dimensional gauge coupling constant and 
the size (and also the shape) of the internal six manifold, respectively. For an isotropic six 
manifold, we have (Re(7})) ^ (Re(T)) where Re(T) denotes the overall modulus whose 
VEV corresponds to the radius- squared of the internal six manifold in the heterotic string 
length unit. The pseudoscalar components Im(S') and Im(7/) are often called the model- 
independent axion and the model-dependent Kahler axions, respectively |[TO| , p!Tll . These 
axion-like moduli are periodic variables and we normalize them by imposing the periodicity 
conditions: 

Im(5) = Im(^) + 1, Im(T/) = Im(T/) + 1. (1) 

Under this normalization, the scalar components are given by ||12[ 

Re(5) ^ ^ ^ 47re-2^-^ ^ 2(4vr«:2)-2/V , 
acuT (2a') 

cl/3 T/l/3 

Re(r) ^ ~ 6i/3(4^^2)-i/3^^^i/3 ^ (2) 
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where e'^^ is the heterotic string dilaton, V is the internal space volume, is the eleven- 
dimensional gravitational coupling constant, and finally np denotes the length of the eleventh 
segment in M-theory. The above relations, together with Mp = 27TpK,~'^V, show that the 
heterotic string coupling e^^ scales as aGt/T[Re(T)]^ and the eleventh length p scales as 
Mp^[Re(T)]3/^ Inserting the proper numerical coefficients, it is easy to see that for acuT ~ 
1/25 and Mp ^ 2 x 10^^ GeV, the large radius limit with (Re(T)) ^ 1 corresponds to 
the M-theory limit with a strong heterotic string coupling, which can be described by an 
eleven- dimensional supergravity on a manifold with boundary HTsj . 

For the superfield normalization determined by the periodicity condition (|l]), the holo- 
morphic gauge kinetic functions satisfy the relation |jl2| 



ATTfE.-ATTfE'^^Y.ll'^I^ (3) 

where //'s are integer coefficients, and /^g and /e'^ denote the gauge kinetic functions for Es 
and i?g, respectively. Here the gauge kinetic functions are normalized as (Re(/a)) = l/gl and 
(Im(/a)) = Oa/Sn'^, where ga and 9a denote the gauge coupling constant and the vacuum 
angle for the a-th gauge group, respectively. Note that for integer // the relation is 
consistent with the periodicity condition and the periodic vacuum angles 9a = Oa + 27r. 
The relation (^) suggests that, even in the M-theory limit, Re(T) ^ Re(7/) can not be 
arbitrarily large, but is constrained not to significantly exceed 47rRe(/£;g) ~ acur ' "^^^ 
compactifications on a smooth Calabi-Yau manifold with vanishing E'^ field-strength, it 
turns out that at least one of //'s is positive and other //'s are still non-negative integers 



T^ . This then leads to the upper limit: 



(Re(r)) < (4) 

acuT 



which corresponds to the lower limit on the Newton's constant Gn discussed in and 
also to the lower limit on the Kaluza-Klein scale Mkk ~ c^gut/VGn discussed in [p!5| . 

The Kahler potential of the effective supergravity model depends upon Re(S') and also 
Re(Tj) with unsuppressed coefficients of order unity. As a result, once the four-dimensional 



supersymmetry is broken, Re(5') and Re(7/) receive the masses of order 7713/2 from the 
supergravity scalar potential. However the masses of the axion-like moduli Im(S') and Im(77) 
are constrained by the approximate non-linear global U{1) symmetries defined as 

U{l)s : Im(5) ^ Im(^) + as, U{l)i : Im(T,) ^ Im(T,) + aj, (5) 

where as and aj denote arbitrary real constants. In the limit where one (combination) 
of these ?7(l)'s becomes an exact symmetry, the corresponding axion-like field becomes an 
exact Goldstone boson and thus is massless. In string/M-theory, these [/(l)-symmetries are 
explicitly broken either by the Yang-Mills axial anomaly or by the world-sheet (membrane) 



instanton effects |jT6[. If the hidden sector gauge interactions provide a dynamical seed for 
supersymmetry breaking, the ?7(l)-breaking by the hidden sector Yang-Mills axial anomaly 
is so strong that the linear combination of Im(S') and Im(7/) which couples to the hidden 
sector anomaly get the masses of order 171^/2. However most of the known compactification 
models allow a combination ofU{l)s and U (1)/ which is free from the hidden sector anomaly, 
but still explicitly broken by the world-sheet instanton effects and/or by the observable sector 
anomaly (mainly the QCD anomaly) [|TT],|TB|. 

Let 0, being a linear combination of Im(S') and Im(7/), denote the axion-like moduli for 
the combination of U{l)s and U{l)i which is free from the hidden sector anomaly. Then 
the effective potential of this axion-like moduli includes first of all the contribution from 
the world-sheet instanton effects which is estimated to be |^|12| 



Vws ^ e-^-<^^(^))mI/2M|. (6) 

If couples to the observable sector QCD anomaly, would include also the contribution 
from the QCD anomaly, Vqcd ~ fl'^'i- ia.ct, one could argue based on supersymmetry 
and the periodicity condition (4) for the axion-like moduli that there is no other type of 
contribution to |T^ . 

For compactifications with (Re(T)) ^ 17, is dominated by the world-sheet instanton 
contribution (|^). (Here and in the following, we assume ^3/2 ~ 1 TeV for the simplicity of 
the discussion.) Then the mass of the axion-like moduli field d) is estimated to be 



^ e-'^^^^^^^^ms/s. (7) 

However if (Re(T)) ^ 17 and also couples to the QCD anomaly, is dominated by the 
QCD contribution, leading to m<^ ~ /^m7r/Mp. (If (Re(T)) is even bigger than about 20, 
we have Vws ~ IO^^Vqcd and then this (j) can be identified as the QCD axion solving the 
strong CP problem PJT^.) In models with hi i > 1, there can be a combination of U {l)s and 



f/(l)/ which is free from both the hidden sector Yang-Mills axial anomaly and the observable 
sector QCD anomaly. The axion-like moduli field for this combination does not couple to 
the QCD anomaly and then its mass is given by Eq. (0) even for (Re(T)) ^ 17. 

The above discussion, particularly Eq. (^, implies that is highly sensitive to the 
compactification radius which is measured by (Re(T)). It can be extremely smaller than 
m3/2 if the compactification radius is large enough to have (Re(T)) ^ 1. As was noted in 
the discussions above (Q), a large value of (Re(T)) is allowed in the M-theory limit but it 
is constrained not to significantly exceed ^^^^ ~ 25. Then for the range < (Re(T)) ^ 25, 
7710 can be anywhere between 7773/2 and the extremely small mass 10~^^m3/2. Note that even 
when (Re(T)) ^ 1 for which the weakly-coupled heterotic string theory provides a sensible 
description for the dynamics above Mp, the axion-like moduli mass can be smaller than 
7773/2 by one or two orders of magnitudes. 



III. MODULI DYNAMICS IN THE EARLY UNIVERSE 

In the previous section, we have noted that, in string/M-theory with a large compacti- 
fication radius, some of the axion-like moduli can be much lighter than the gravitino mass. 
Motivated by this observation, in this and next sections, we study cosmological aspects of a 
generic light modulus with an arbitrary mass below the weak scale. 

Moduli dynamics in the early universe would be governed by the free energy density V^g 
which depends not only on (p but also on other scalar fields $ and the radiation temperature 
T. Expanding Ves around the present moduli VEV which is set to zero, one generically has 
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Kff(0, T) = fio($, T) + T)0 + l^]2($, T)02 + (8) 

where the moduh-tadpole Qi arises since the expansion is made around the present moduh 
VEV, not around the minimum of Ves in the early universe. Obviously at present, f2„ are 
given by 

noi%, To) = (]i($o, To) = , (]2($o, To) = ml , (9) 

where $0; To, and denote the present VEV of the present temperature, and the 
present mass of 0, respectively, and the vanishing of f2o('^'o;To) comes from the vanishing 
(or negligibly small) cosmo logical constant at present. 

In the early universe, $ and T can take values far away from the present ones. As 
a result, fi„'s in the early universe can significantly differ from their present values. For 
instance, when H ^ m^, the effective moduli mass y/Tl^ can be of order H and thus much 
bigger than m^. 

To proceed, let us parameterize the free energy density as follows P^ : 



Kff(0, T) = + \^H\ct^ - 00' + ... , (10) 

where denotes the moduli mass at present and the ellipsis denote the irrelevant terms. 
This parameterization is useful since in most cases of interest, c and 0i are approximately 
time-independent constants, which would greatly simplify the analysis. Obviously the coef- 
ficients c and 01 describes the non-derivative interactions of with the background energy 
density in the early universe, e.g. the radiation or the infiaton energy density. Generically 
01 is expected to be of order Mp. The value of c measures how strongly couples to 
supersymmetry breaking environment and thus 

(11) 

It was noted in Ref. [0 that in the case of c ^ 1, the initial moduli misalignment can be 
rapidly damped away. In this paper, we consider only the case c ^ 0{1) which appears to be 
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more natural. Note that for axion-like moduli discussed in the previous section, c ~ g-i-ReC^) 
can be very small. 

Treating c and (pi as constants, the equation of motion for is given by 

+ 3H<j) + {ml + c'^H^)(f) = c^H^(f)i . (12) 

In the following, we wish to examine the evolution of by solving this equation with the 
initial value taken at some initial time tf. 

0(t,) = 4>i , <P{U) = . (13) 

Let us first consider the evolution of during inflationary era for which the Hubble 
parameter H is roughly a constant. The solution of Eqs. (^2|) and ( pTSD is easily found to be 
as follows: For (f-ff)^ > m| + c^H'^, 



4>{t) = 0mm + (0i - 01 

while for (|i/)^ < "^i + c^H'^, 



2f3 2(3 



(14) 



0(t) = 0min + (0i - 0min)e-^^(*-*°) 



3H 

cos[p'm^{t - to)] - — sm[p'm^{t - to)] 



, (15) 



where P = yl — |(c^ + j^), P' = y'l — (| — c^)^-, and 0mm is the temporal minimum of 
the effective potential which is given by 

0mm = • (16) 

Let us consider some interesting limits of the above solution. If <C H and c -C 1, we 
have /5 ~ 1 — |(c^ + -^) and then the modulus value after the inflation (0/) is given by 

(0.-0.m)(^^^jf^), (17) 

where is the number of e-folding. In the above, we assumed that c and m^/H are small 
enough so that A^'e -C ^37^7^- Thus simply speaking, in the case of <^ H and c <^ 1, 
is frozen at its initial value 0^. This is what we thought to happen during inflation before 
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realizing the role of H"^ term in Eq. (|T0|). But for <^ H and c ~ 1, we have the final 
modulus value 

0/ ~ 0min + (0. - 0min) X Ole-^^^-^/^) , (18) 

showing that rapidly approaches to the temporal minimum ip-cmn — 0i independently of 
the value (pi before the inflation. Finally for ^ H, (p exponentially approaches to 
0min — ^^^01, but with an exponentially decreasing oscillatory tail. This is what happens 
in so-called thermal inflation p| which will be discussed in more detail in the next section. 

Let us now consider the evolution of during the radiation-dominated (RD) or matter- 
dominated (MD) era for which the Hubble parameter is given hj H = p/t where p = ^ 
(RD), I (MD). In RD or MD era, the solution of Eq. (|l|) is given by 

0(.) = p\^<f>,^^^^ + Ci^ + , (19) 

where z/^ = — p'^c^ > for a = |(3p — 1) and z = m^t. Here J^{z) and Yi,{z) are Bessel 
functions, and S^^y{z) is Lommel function which is defined by 



Y,{z) r y^Uy)dy-Uz) fy^^Y^dy 
Jo Jo 

/i — ZZ + l /i + ZZ + l f /i — /i — 



+2^''r(^^^ )r(^^^ ) |sin(^7r) J,(z) - cos(^7r)n(^)| . (20) 



The coefficients Ci and C2 are fixed by the initial condition (|T3D: 



C2 = A2cj)i + B2C^i (21) 



where 



A^ = '^z'l{z,Yl{zi)~aY,{zi)} 
M = --zt {ziJlizi) - aJ^{zi)} 

Bl = -^Zi \Yl{Zi)Sa-l,u{Zi) - Y^{Zi)S'^_;i^ „{Zi)^p 
71 
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B2 = -z, {jUz,)S^.,,,{z,) - Mz,)S'^_,^,{z,)}p'. (22) 



2 

10 



Here Zi = rUfpti and the prime denotes the differentiation with respect to z. 



From the solution (0), we can calculate the moduli abundance coming from the coherent 
oscillation. For z ^ 1, (j){z) is dominated by the oscillating tail 

m ^ {lY' ^-^^ {Ci cos[z - (z. + 1)^] + C2 sin[z - + ^)^]} • (23) 

The energy density of this oscillating modulus is given by 

P, = + 0^) - [I) {Cl + Cl)mlz-'^ (24) 

and normalizing it by the entropy density s = ^g^T^, we find 

^<^n^'-^ = ^iC! + Cl)^^, (25) 

where g^, denotes the effective number of the relativistic degrees of freedom at T. 

Most of the cosmological implications of (j) is in fact determined by the oscillation am- 
plitude: 

(ci + ci 



TT 



(26) 



which we call the initial moduli misalignment throughout this paper. The oscillation coeffi- 
cients Cl and C2 are determined by the modulus value 0, = (j){ti) at an initial time ti with 
(/)'(tj) = 0, and also by the two dynamical parameters c and 0i in the free energy density 
(p!0|) which governs the moduli dynamics at later time t ^ ti. Depending upon the periods 
under consideration, Zi = m^ti may be chosen to be either very small, or of order unity, or 
very large. It turns out that the coefficients Ai^2 and -81^2 in Eq- (|2^ ) are essentially of order 
one for -Zj ~ 1 and an arbitrary value of c ^ 1. Thus roughly speaking, Ci and C2 are linear 
combinations of (pi (at ti ^ m'^'^) and c^0i with coefficients of order one. At any rate, using 
Eqs. ( PUD and (^), any cosmological bound on S(f) can be translated into a constraint on the 
parameter set (0j,c^0i). 

If the universe were radiation-dominated when 0-oscillation begins at t ~ m^^ and there 
is no entropy production since then, we have z = (^ 2n^g ) ~^t^^ and thus 
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'"^^■*-(2^)'73'""*''^'"^- '''' 

In other case that the universe was matter-dominated, for instance by the inflaton oscillation 
or by the heavy moduli oscillation, at the moment when 0-oscillation begins, one has to take 
into account the subsequent entropy production due to the out-of-equilibrium decays of the 
inflaton or heavy moduli. Assuming that the whole matter energy is converted into the 
radiation with the reheat temperature Tr, the modulus energy density normalized by the 
entropy density is given by Eq. (|25|) evaluated at the reheat time tn = ^ 



IV. THE CONSTRAINTS ON THE INITIAL MISALIGNMENT 

A. The cosmological moduli problem 

In the previous section, we obtained the moduli energy density normalized by the entropy 
density when the coherent oscillation begins during the RD era: 

. . a X 10= (J^ f )^ GeV . (29) 

A coherently oscillating modulus can dominate the energy density of the universe unless 
m^Yff, is less than the temperature T^q ~ 3eV of matter-radiation equality. This implies 
that one has to worry about the over-produced moduli which would contradict with the 
cosmological observations as long as the moduli mass is in the range 

m<^> 10-26 eV. (30) 

On the other hand, sufficiently heavy moduli decaying before about one second do not affect 
the standard prediction of the big-bang nucleosynthesis, and thus would not contradict with 
the currently known cosmological observations. Having interactions suppressed by Mp, the 
moduli lifetime is estimated to be 
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-0-exlO-(^)" sec, (31) 

where ^, being roughly of order one, is a coefficient which accounts for the ambiguity in 
our estimate of the hfetime. For moduh hghter than ~ 40 TeV, their decay products may 
change the abundance of the hght elements after or during the nucleosynthesis, the spectrum 



of cosmic background radiation, or the observed 7 and X-ray backgrounds []T9|. Moduli 
whose lifetime is longer than the age of the universe would overdose the universe unless 
m^Y^ ^ 3eV. This means for the moduli mass ^0.1 GeV, the initial misalignment is 
constrained as 

unless there is a late entropy production after cj) starts to oscillate. 

Based on the analysis of |1^, the recently reported X-ray background ||20| and Eq. (p2D, 



we obtain the constraints on 50 arising from these considerations again under the assumption 
that there is no entropy production after 0-oscillation begins at t ~ ■ The results 
are summarized in Figure 1 showing the cosmological upper limit on the initial moduli 
misalignment 5(f) for the moduli mass below ~ 40 TeV. The line (0) comes from Eq. (|32D 
which is required for the moduli not to overdose the universe, and the line (1) from the 
recently reported X-ray background. The lines (a)-(h) are obtained from the observed 7-ray 
background, the spectrum of cosmic microwave background radiation, and also the light 
element abundances. 

The results of Figure 1 show that, in the absence of a late time entropy production, the 
initial moduli misalignment is required to be very small compared to its natural value ~ Mp, 
for instance 50 < 10~*^Mp for ^ 1 eV, 50 < lO^^^Mp for ^ 1 MeV, 50 < lO^^Mp 
for 1 GeV. Let us recall that 50 is determined by the modulus value 0j at an initial 

time ti with 0'(tj) = 0, and also the dynamical parameters (c, 0i) in the free energy density 
which would govern the moduli dynamics at time t ^ tj. (See Eqs. (^T])-(^B]).) Its order 
of magnitude is roughly given by the bigger one among 0j and c^0i for the initial time 
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ti ^ Tn^^. Its natural value would be of order Mp which then falls at far above the upper 
limits in Figure 1 for most range of the moduli mass. 

B. Dilution by heavy moduli decays 

If there is an entropy production during the period after (f) starts to oscillate but well 
before decays, the moduli energy density is diluted as — * '^4>/^ where A = Rafter /^before 
denotes the entropy production factor. Since oc 50^, this obviously leads to the relaxation 
of the constraints on the initial moduli misalignment which will be discussed below. 

Usually the entropy production is due to out-of-equilibrium decays of non-relativistic 
particles which appear in the form of another coherently oscillating scalar field ip. In order 
to be compatible with the big-bang nucleosynthesis, this entropy-producing scalar field if is 
required to decay before the nucleosynthesis with the reheat temperature Tr > 6 MeV pT 



Also note that if the light moduli is so light that its oscillation begins after the entropy 
production by (f is over, i.e. ^ gl^'^T^/Mp, is not affected by the entropy production 
by ip. Thus the energy density of oscillating can be diluted only for the moduli mass 

To be more specific, let us consider the interesting possibility that the entropy-producing 
field <y9 is a massive moduli (f with 6(f ~ Mp and ~ 40 TeV which would give a maximal 
entropy production with Tp ^ 6 MeV. This heavy modulus dominates the energy density 
of the universe as soon as its coherent oscillation begins at T ^ 10^^ GeV. The light moduli 
which are lighter than ~ 40 TeV but heavier than the bound in Eq. (^) start to coherently 
oscillate during the matter dominated era by the heavy modulus oscillation. Then we can 
apply Eq. (pSD and get the abundance 

which is diluted compared to Eq. ( P7D by a factor Rafter /5'before = {fnipMpY^'^ / gl^^Tp. It is 
rather straightforward to derive the cosmological limits on the light moduli misalignment 
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(50 for the diluted moduli density (|3^) as we did for the case without any late entropy 
production. The results for the case of = 6 MeV, i.e. maximal dilution, are shown in 
Figure 2. It shows that the entropy production by heavy moduli decays ameliorates but 
does not completely solve the cosmological moduli problem. 

C. Dilution by thermal inflation 

The most efficient way to dilute the dangerous light moduli (f) is to have a late inflation 
since inflation dilutes the moduli density both by the spatial expansion and by the large 
amount entropy production. The most natural framework for a late inflation would be the 
so-called thermal inflation models P]. 

In thermal inflation models, the entropy-producing fleld (p corresponds to a flaton fleld 
parameterizing a flat direction in supersymmetric models. This flat direction is lifted by the 
soft breaking mass and also by the Planck scale suppressed non-renormalizable terms, leading 
to the following (renormalization group improved) effective potential at zero temperature: 

K, = K,-mJ|^r + ^^, (35) 

where n is a model-dependent integer and the negative mass-squared can arise as a conse- 
quence of radiative corrections associated with the strong Yukawa coupling of (p. The true 
vacuum expectation value is (y?) ~ (m^Mp)^/'-"^"'^^ and Vq is adjusted to Vq ~ m^(y9)^ ^ 

2n+4 2n 

m^^^ Mp^^ in order for the true vacuum energy density to vanish. At high temperature 
T ^ m^, the effective flaton mass-squared including the thermal contribution of (9(T^) is 
positive and thus {^)T-»m^ = 0. For the period of ^ T < Vq^'^ ^ the universe is vacuum- 
dominated, yielding an exponential expansion with the number of e-foldings (A^^e) which is 
determined by 6^"= ^ 0.42(100/5(*)i/^(Mp/m^)"/2"+2_ At T ^ m^, the flaton starts to roll 
down to its true vacuum value and then the vacuum energy density Vq is converted into 
the energy density of coherently oscillating flaton fleld. The oscillating flatons eventually 
decay and are converted into the radiation with the reheat temperature 
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n-1 



Tn - 1.7^-/^ - 0.l7^/^^;^^/V (^) , (36) 

where 7 is introduced to parameterize the yj-decay width = 'ym^/64:TT{ip)'^ and denotes 
the effective number of the relativistic degrees of freedom at Tr. As was discussed in cp 
can couple to ordinary matter through its mixing with Higgses, and then the most efficient 
decay channel is the decay into stop and anti-stop pair. Assuming > 2mj so that (p can 
decay into stop pairs t and t*, we have roughly 7 ^ (2mf/m^)^. Although it can be a quite 
small number (particularly when ^ m^), we assume here 7 ~ 1 as a conservative choice. 
The entropy production factor of thermal inflation is given by 



O.ll-i ( ^ ) . (37) 



2n + 2 



S'before ^TrITI^ ' / 

From Eq. ( P^D we see that to achieve Tji ^ 6 MeV it is required that 7^/^m<^ > 100 MeV, 
7^/'^m<^ > 60 GeV, and 'y'^/^m^ ^ 700 GeV forn = 1, 2, 3 respectively. Under this restriction, 
the maximal entropy production factor is 2 X 10l^ 6 X 10^371/7^ 8 X 102^7^/^ respectively. 

The light moduli can start to oscillate either before or after the thermal inflation 
depending upon their masses. Obviously the moduli oscillation should start before thermal 
inflation occurs in order for thermal inflation to sufficiently dilute the moduli density. This 
gives a lower bound on the moduli mass for which thermal inflation leads to the sufficient 
dilution of moduli density: > Hti ^ {Vo/M],y/^ ^ 100 eV, 100 keV, 1 MeV for n = 
1, 2, 3 respectively, where Hn is the Hubble expansion parameter during thermal inflation. 
The moduli lighter than this bound but start to oscillate before the flaton decay are also 
diluted somewhat by the entropy production due to the flaton decay. For this mass range the 
analysis of the previous section can be applied. For ^ iTxi? we need a detailed analysis, 
but we will not concern such details. In Figure 3(a), we showed the relaxed constraints 
on the initial misalignment 6(f) when the moduli density is maximally diluted by thermal 
inflation. 

Although thermal inflation provides a huge entropy and thus dilute the moduli energy 
density due to the intial misalignment, it can cause an additional misalignment induced 
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by a shifted minimum of the free energy density during the inflation. As was noted in the 
discussions below Eq. (|18|), during the thermal inflation period, the shifted minimum is given 
by 0min ~ —p-4>i- Foi' -^Ti ~ "^(^5 this shifted minimum leads to an energy density 

where c and 0i are those in the free energy density ([T0| ) during thermal inflation. We 
showed the resulting constraints on 0i in Figure 3(b) for the case that c ~ 1 over the entire 
moduli mass range, and in Figure 3(c) for the more plausible case that c ~ '^<^/'^3/2 with 
7713/2 = 100 GeV. 

V. CONCLUSION 

In this paper, we have discussed the possibility of light moduli having a mass far below 
the weak scale, and examined the cosmological bounds on the initial moduli misalignment 
for the mass range 40TeV ~ 10^^^ eV. A very light moduli can arise as an axion-like 
moduli in string/M-theory with a large compactification radius, ^ e"'^^''*^^-* 7773/2 with 
Re(T) = 1 ~ ^^^^ . Also generic moduh in gauge-mediated supersymmetry breaking models, 
can have a mass in the range 10 eV ~ 1 GeV. We then studied the cosmological evolution 
of a generic light modulus to quantify its relic energy density which is determined by the 
initial misalignment The initial misalignment 5(f) is set by the modulus value (pi at an 
initial time tj with (j)'{ti) = and also by the dynamical parameter c^0i in the moduli free 
energy density ([10|) at t > t^. 

For the case that there is no entropy production after 0-oscillation begins at t ~ 777^"^ 
(but of course before decays), the bounds on the misalignment 6(f) coming from various 
astronomical and cosmological observations are shown in Figure 1. We then discussed how 
much such bounds can be relaxed by a late entropy dumping which can be driven typically 
by e.g. out-of-equilibrium decay of another heavy moduli [Figure 2] or by a late vacuum 
domination like thermal inflation [Figure 3]. The bound on the initial misalignment turned 
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out to be severe so that 54) ^ Mp is not allowed in most cases. Thermal inflation appears 
to be efficiently relax the bound, however still Scf) fa Mp is allowed only for a limited moduli 
mass range as shown in Figure 3. 
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FIGURES 

FIG. 1. Constraints on the initial misalignment of moduli. We take = ^, h = 0.7 and the 

lifetime r = 10^^ ( iGeV ) The line (a) comes from the lower bound on the age of the universe, 

^(l,h^ < 1, which is unavoidable even for very light moduli. The lines from (b) to (j) are for the 
case that moduli decay to photons or charged particles: (b) Recently reported X-ray background, 
(c) Observed 7-ray background, (d) Distorsion of CMBR, fi > 8 x 10^"^, (c) Photoproduction, 
(D+3He)/H > 10"^, (f) Photodestruction, D/H < 10-^ (g) Hadronic shower, (D+3He)/H > 10"^, 
(h) Hadronic shower, i^(^He) > 0.25, (i) Entropy production, (D+^He)/H > 10~^, (j) Entropy pro- 
duction, YpC^He) > 0.25. The line (a) can be straightforwardly extended to the region < 1 eV. 
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FIG. 2. Relaxed bounds on the initial misalignment of light moduli when there is an entropy 
production due to the heavy modulus decay which yields the reheating temperature Trh = 6 MeV. 
Again the line at ~ 1 eV can be straightforwardly extended to the region < 1 eV. 
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FIG. 3. Constraint on the misalignment of moduli when there exists thermal inflation with 
Tr = lOMeV: (a) Constraint on the initial misalignment of moduli, (b) Constraint on the mis- 
alignment {S(f) = <f>i) induced by the shifted minimum during thermal inflation when c = 1 over 
the entire mass range, (c) Constraint on the misalignment {d(f) = ^i) induced by thermal inflation 
for c = with = 100 GeV. 
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